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ABSTRACT 

Measurements  of  the  ratio  of  the  deuteron- to- proton  electric  form 
factors,  G_  /G_,  ,  were  made  from  elastic  electron- deuteron  scattering 

Hi ,  iii 

d  p 

to  a  precision  of  approximately  1$  for  the  range  of  momentum  transfers, 

O  rJ  r\ 

0.10  S  q  ^0.8  f“‘~  ,  and  for  electron  scattering  angles  of  45°  to  120°. 

It  was  found  that  within  experimental  errors  the  slope  as  obtained  from 
the  ratio  G_,  /G„  agrees  with  the  extrapolated  thermal  neutron-electron 

XL,  ill 

n  p 

interaction  slope  when  relativistic  corrections  and  Feshbach-Loraon 
deuteron  wave  functions  are  applied  to  the  electron-deuteron  scattering 
results.  The  deuteron  radius  was  found  to  be  1  ..95  -  0.02  f  ,  which  agrees 
with  the  radius  predicted  by  the  Feshbach-Lomon  calculation. 

*  *  * 

We  have  measured  the  ratio  of  elastic  electron-deuteron  scattering 

to  elastic  electron-proton  scattering  in  the  range  of  momentum  transfers 
o  .2 

0.1  S  t1"  •  0.8  f  with  a  precision  of  one  percent  or  better. 

It  is  possible  to  extract  from  such  measurements  values  for  the 
ratio  G„  /G„  .  In  our  range  of  momentum  transfers  this  rati...  is  very 

Jci  h, 

n  p 

close  to  the  value  of  G_  .  The  proton  charge  form  factor  G„  has  r.ot 

E  r  E 

n  p 


1 

.L 


been  extensively  measured  at  such  low  values  of  q  and  the  use  of 
experimental  fits  to  data  at  higher  momentum  transfers  to  extract  G- 

I It 

n 

must  be  made  with  the  assumption  that  there  are  no  unusual  fluctuations 
in  this  region.  We  report  the  ratio  G_,  /G  and  apply  the  best  known 

bt  I* 

n  p 

fits  of  Gg  to  extract  a  value  of  G£ 

P  n 

Few,  but  precise,  data  have  been  reported  in  our  range  of  momentum 
transfers.^-  This  work  has  been  done  to  look  more  accurately  at  the 

discrepancy  between  the  results  of  neutron-electron  Interaction  measure- 

2  2  2 
ments,  which  yield  a  positive  slope  dG_  /dq  at  q  -  0  ,  and  the 

n 

fact  that  the  value  and  the  slope  of  G..  as  extracted  from  electron- 

L 

n  2 

deuteron  scattering  averaged  to  zero  at  all  values  of  q  .  This  apparent 
disagreement  was  not  removed  by  previous  measurements  at  low  momentum 
transfers. ^ 

The  extraction  of  G£  from  electron-deuteron  scattering  data 

n 

requires  knowledge  about  the  deuteron  wave  functions.  Major  efforts  have 
been  made  in  the  last  ten  years  (see,  e.g..  Refs.  3,  4,  and  5)  to  improve 
the  status  of  the  theory  of  the  deuteron.  Gross^  has  shown  that  relati¬ 
vistic  corrections  to  the  wave  functions,  although  small  in  the  range  of 
2 

q  considered  here,  can  play  a  very  important  role  in  the  extraction  of 

G_  .  Casper  and  Gross ^  found  that  with  the  proper  choice  of  realistic 
n 

wave  functions  and  the  application  of  relativistic  corrections  it  was 
possible  to  affect  the  data  in  a  direction  that  would  tend  to  reduce  the 
difference  between  the  neutron-electron  interaction  results  and  those  from 
electron-deuteron  scattering. 

The  experiments  were  performed  at  the  linear  accelerator  of  the 

2  -2 

Naval  Postgraduate  School  at  Monterey  in  the  range  0.1  $  q  i  0.4  f  and 


2 


at  the  Mark  III  accelerator  at  Stanford  University  in  the  range 
2  *2 

0.?  S  q  r  OiJ  f~r  .  The  experimental  cet-up  was  basically  the  same  in 
both  laboratories. 

A  well-defined  beam  of  electrons  with  a  momentum  spread  of  §  0,1?S 
passed  through  the  targets  and  then  into  a  beam  collector.  The  targets 
used  were  polyethylene  (CH^Jn  and  deuterated  polyethylene  (CD^Jn  . 

A  carbon  target  was  used  for  background  subtraction.  The  target  thick- 
nesses  (<  10  J  rad.  lengths)  were  matched  for  equal  energy  loss  by 
ionization  to  optimize  the  conditions  for  the  carbon  subtractions.  The 
normalization  factor  for  this  subtraction  was  obtained  by  measuring  the 
number  of  counts  under  the  carbon  peaks  in  the  three  targets  at  the 
various  energies  and  angles.  The  scattered  electrons  were  analyzed  in 
momentum  in  a  magnetic  spectrometer  and  counted  in  a  multichannel  array 
of  scintillator  counters.  The  relative  efficiencies  of  the  counters 
were  separately  determined  before  each  set  of  data  points. 

The  absolute  energy  of  the  incoming  electrons  was  known  to  «  0.2$ 
and  monitored  by  nuclear  magnetic  resonance.  The  spectrometers  were 
monitored  and  stabilized  by  accurate  rotating  coil  fluxmeters. 

The  scattering  angle  at  the  Stanford  laboratory  has  been  determined 
to  0.05°  and  to  0.1°  at  the  Monterey  laboratory.  It  is  clear  that  many 
parameters  that  would  influence  an  absolute  measurement  cancel  out  in  a 
ratio  determination,  e.g.,  solid  angle,  charge  integration,  absolute 
counter  efficiency,  etc.  To  this  end  we  have  kept  the  scattering  angle, 
solid  angle  and  scattered  momentum  fixed  for  the  carbon,  CH„  and  CD0 

c  d 

measurements  at  each  data  point.  The  incident  energy  was  adjusted  to 

2 

give  the  correct  q  .  Only  small  changes  in  energy  were  necessary,  as 


shown  in  fable  I. 


Radiative  corrections  were  calculated  according  to  Tsai. '  Because 
of  a  possible  contribution  by  quasi-elastically  scattered  electrons  from 
t,he  deuteron,  cut-offs  on  the  radiative  tails  of  the  deuteron  and  the 
proton  peaks  were  made  at  energies  less  than  2.2  MeV  below  the  peak. 

We  have  shown  that  the  ratios  formed  from  our  data  are  independent  of  the 
cut-off  position  in  the  radiative  tail. 

The  deuteron  form  factor  is  determined  experimentally  by 


\dft/ 


Exp 


,  clo-\ 
vdft/ 


X  G 


Mott 


D 


(1) 


2  6 
Gp  can  he  brought  into  the  form 


gH  =  A(q?)  +  B(q2)  tan2  e/2  . 


(2) 


With  q  = 


q 


4m' 


«  1  in  our  range  and  making  1  +  t|  ->  1  ,  we  have 


A(q2)  =  G2  +  |  rf  G2  +  |  q  G2  and  B(q2)  =  3  q  G2  .  (3) 

G_  ,  G_  ,  and  G,,  stand  for  the  physical  charge,  quadrupole,  and  magnetic 
form  factors  of  the  deuteron,  respectively.  They  are  connected  to  the 
physical  nucleon- elec trie  form  factors  for  the  proton  and  neutron  (G_ 


and  G^  ,  respectively)  by 

n 


GC  ■  (0E  *  °E  1  rc 

P  n 


(**) 


°«  -  <°E  *  aE  >  D« 

p  n 

Gm  =  (G„  +  G.,  )  +  (GM  -  GM  )  d[J 

M  E  E  M  M  M  M 

p  n  p  n 


E  M 

The  Dc  ,  D,  ,  D.„  ,  and  are  integrals  over  suitably  chosen  wave- 


4 


functions  of  the  deuteron.  The  model  dependence  of  our  result  is  contained 
in  these  terms. 

In  the  kinematic  range  of  this  experiment,  =  p^E^  is  a  good 

a 

approximation.  The  expression  for  G^,  is  considerably  simplified  when 

this  is  assumed,  as  well  as  "scaling"  for  the  proton  and  neutron;  i.e., 

G„,  =  G„  and  G.,  =  (j  G„  .  (The  p's  axe  the  magnetic  moments  of 
M  p  E  M  n  E  ° 

P  P  n  p 

the  respective  particles.) 

The  results  of  the  Drickey  and  Hand'*'  experiment  directly  support 

this  assumption  in  the  momentum  transfer  range  of  our  experiment.  They 

find  Gg  =  G^  /p  =  1.008  ±  .014  .  In  the  deuteron,  our  maximum  magnetic 
P  P  P 

correction  is  <  1$  ,  so  that  a  1#  violation  of  scaling  has  a  negligible 
effect  on  our  results. 


1 


(8) 


Dc  GE 


In  our  range  of  momentum  transfers,  the  quantities  D,  and  are 

P 

close  to  one,  so  that  a  small  change  in  either  of  them  has  a  strong 

influence  on  G„ 

E 

n 

Much  of  the  earlier  experimental  work  on  elastic  electron-deuteron 

Q 

scattering  has  been  analyzed  with  wave  functions  developed  by  Partovi 
from  a  Hamada  potential  with  a  core  radius  of  0.48-5  F  and  a  D- state 
contribution  of  7$* 

Feshbach  and  Lomon^  proposed  a  deuteron  boundary- condition  model 
with  a  core  radius  of  0.735  F  and  a  D- state  value  equal  to  4.9$. 
(Within  the  limits  from  4.6$  and  6.1$  their  values  of  vary  only  by 

about  1  part  in  10  J  in  our  range  of  momentum  transfer.) 

Many  attempts  have  been  made  to  derive  a  relativistic  theory  for 
elastic  electron-deuteron  scattering.  Grossb  made  the  most  complete 
attempt  to  find  relativistic  effects  by  treating  the  deuteron  wave 
function  in  a  relativistic  manner.  These  calculations  describe  the 
distortion  of  the  wave  function  of  a  moving  deuteron.  Corrections  stem¬ 
ming  from  relativistic  modifications  of  the  nucleon  current  seem  to  be 

7 

small  and  are  neglected.  Casper  and  Gross  derived  a  relativistic 
treatment  of  the  wave  function  that  yields  an  additive  correction, 

/v.  ,  to  G„  .  Such  a  correction  is  approximately  equal  to  q  /8m~ 

P 

n  n 

and  it  is  reasonably  model  independent.  They  also  showed  that  meson 

2  -2 

exchange  effects  give  negligible  contributions  below  q  «  10  f 

The  Partovi  (P)  and  the  Feshbach- 1  omor.  (FL)  models  with  and  without 
relativistic  corrections  are  applied  to  our  data  in  the  next  section. 


6 


In  Table  T  we  present  our  experimental  data.  The  error  quoted  con¬ 
sists  of  counting  errors  (including  the  carbon  subtraction  and  efficiency 
corrections)  and  an  error  of  in  the  target  thickness  determination. 

Table  II  shows  four  sets  of'  values  of  G_  /G„  as  obtained  with 

n  p 

(a)  the  Feshbach-Lomon  wave  function  with  relativistic  corrections, 

(b)  the  Feshbach-Lomon  wave  function  without  relativistic  corrections, 

(c)  the  Partovi  wave  function  with  relativistic  corrections,  and 

( d)  the  Partovi  wave  function  without  relativistic  corrections. 

Values  for  the  relativistic  correction  £G  and  the  structure 

hi 

n 

factor  Dc  for  the  Feshbach-Lomon  wave  function  are  also  given. 

p 

In  each  of  the  columns  we  show  the  slope  d(3g  /  G£  )/dq  to 

n  p 

illustrate  the  rapid  decrease  in  the  slope  when  going  from  a  reiativistically 
corrected  Feshbach-Lomon  wave  function  to  an  uncorrected  Partovi  wave 
function.  These  fits  assumed  a  value  of  zero  for  the  intercept.  (We 
are  aware  of  the  slight  inconsistancy  in  adding  to  the  ratio 

ge  /°e  • ) 

n  p 

In  order  to  extract  values  of  G_  we  have  used  the  deVries  b' 

E 

11  n 
fit  for  the  proton  form  factor  G„  .  This  fit  is  in  very  good  agree- 

EP 

ment  with  the  absolute  measurements  of  G„  that  were  made  by  Driekey 

p 

arid  Hand.  Table  III  shows  the  results  of  this  application  for  the  two 
extreme  cases  presented  in  Table  II.  The  slopes  as  determined  by  the 
two  evaluations  are  also  given  for  comparison  with  the  neutron-elect ron 
interaction  slope.  The  errors  are  weighted  by  the  X"  of  the  fit. 

It  is  interesting  to  note  that  if  a  non- zero  intercept  were  assumed, 
then  our  results  still  yield  a  consistent  value  of  the  slope.  In  such, 
a  case  we  find  the  slope  =  0.022  ±  0.00b  and  the  intercept  is 

1 


-  O.OOF  t  (  .  k  3  for  the  data  evaluated  with  the  Feshbach-Lomon  wave 

functior  and  the  relativistic  correction. 

We  can  use  oui*  measurements  of  GT,  /G„  to  extract  a  value  of  the 

Eg  E 
d  p 

deut.eroi:  radius,  r  .  However,  an  expansion  of  the  structure  function, 

o  _  n 

,  becomes  model  dependent  for  values  of  q‘  as  low  as  0.1  f  .  In 


the  expansion 


D 

c 


> 
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k 

+  Aq 


(9) 


the  value  of 

.1  £  q?  £  .3 


>.  from  the  Feshbach-Lomon  calculation  is  0.3  in  the  range 

*'  Expanding  Eq.  7  and  including  Eq.  9>  we  have 

o 

(<rd>-<rn>)*X«1‘-  (10) 

d  p 


<  r  >  =  0.116  f  and  is  obtained  from  the  neutron-electron  slope, 

n  ^ 

o  c 

The  best  fit  to  our  data  points  to  Eq.  10  up  to  q^  =  0.4  f"^ 

yields  r  =  1  .°c  f  .OF  f  .  A  change  of  X  by  ±.0n  changes  the 

deuteron' s  radius  by  one  standard  deviation.  A  good  fit  to  our  data 
4 

requires  the  q  t*  rm.  Higher  order  terms  in  might  be  important 

beyond  this  value  of  q  .  The  Feshbach-iomcn  calculation  gives 
r  =  i  f  . 

u. 

In  Fig.  1  we  present  the  contents  of  Tab III.  The  three  points 
of  Drickey  and  Hand  (as  re-evaluated  by  Buchanan1'')  that  fall  in  our 
range  of  q  are  also  plotted. 

Barring  unexpected  fluctuations  of  Gg  in  our  range  of  momentum 

P 

t  ra.usfers,  wr  cor. elude: 


(a)  When  *  hr  Feshbach-Lomon  wave  function  together  with  relativistic 
corivci  Lons  Is  applied  to  the  data,  one  finds  agreement  between  :  he 
neu  .so:  .-e  . nc Iron  slope  at  q  -r  0  and  the  slope  given  by  values  of 
if.  t  he  range  0 .10  :  q  •  0.H0  . 


(b)  Even  within  the  relatively  large  errors  propagated  into  Gp  ,  the 

•  n 

Partovi  wave  function  with  relativistic  corrections  yields  a  slope  that 
is  in  disagreement  with  the  neutron-electron  interaction  slope. 

(c)  There  is  agreement  within  errors  between  the  deuteron  radius  from 
experiment  and  that  obtained  from  the  Feshbach-Lomon  wave  function  that 

was  used  in  (a)  to  find  the  slope  of  G„ 

h 

n 
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table:  captions 


Table  I  Summary  of  experimental  data.  Columns  1,  2,  3>  and  4  contain 

the  momentum  transfer,  scattering  angle  and  incident  electron 

energies  for  the  proton  and  deuteron  data  points,  respectively. 

ColumrPt>  and  6  are  explained  in  Note  b.  Column  7  has  the 

measured  ratios  of  deuteron  electric  form  factor  to  proton 

electric  form  factor.  Column  8  is  the  total  experimental 

error  for  each  point  and  Column  9  is  the  weighted  average 

plus  errors  for  the  points  at  each  value  of  momentum  transfer. 

Table  II  Column  1  is  the  momentum  transfer.  Column  2  is  the  relativistic 

FL 

correction  to  G^.  as  calculated  by  Gross.  Dc  in  Column  3 
n 

is  the  deuteron  structure  function  as  calculated  by  Feshbach- 

Lomon.  The  succeeding  columns  are  the  results  for  G^  /Gg 

n  p 

using  relativistic  corrections  or  not  and  Feshbach-Lomon 
(FI,)  or  Partovi  (P)  structure  functions  for  the  deuteron. 

The  boxes  at  the  bottom  give  the  slope  of  the'  best  line  fit 
to  G_  /G_  for  the  different  cases. 

hi  hi 

n  p 

Table  III  The  de  Vries  b'  fit  is  used  to  find  G„  for  three  of  the 

E 

n 

cases  shown  in  Table  II.  Column  2  gives  the  value  of  Gg 

P 

used. 
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Table  III 


EXTRACTED  SLOPES  dGp  /dq^ 

'n 


2 

q 

(f-2) 

r  a 
ge 

P 

deVries 
b'  fit 

ge 

n 

FL 

°E  +£CE 

n  n 

FL 

ge 

n 

P 

0.10 

0.9881 

-0.0026 

-0.0020  ±  0.0054 

-0.0058  ±  0.0054 

0.20 

0.9765 

40.0039 

+0.0051  ±  0.0031 

40.0017  ±  0.0031 

0.30 

0.9652 

-0.0006 

+0.0012  ±  0.0035 

-0.0037  ±  0.0035 

0.40 

0.9540 

+0.0077 

40.0101  ±  0.0036 

40.0036  i  0.0036 

0.50 

0.9^31 

40.0037 

40.0067  ±  0.0074 

-0.0015  ~  0.0074 

0.60 

0.9321 

40.0011 

40.004 !  ±  0.0066 

-0.0042  ±  0.0066 

0.80 

0.9115 

40.0136 

+0.0164  ±  0.0059 

+0.007 6  -  0.0059 

Straight  line  slope 

(  with 

zero  intercept): 

40.0179  ±  0.0036 

40.00:6  ±  o.ooj.6 

Slope  n  -  e  interaction:  0.0195  -  0.0004 


a: 


No  error  has  been  applied  to 
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FIGURE  CAPTION 


Figure  1 


Straight  line  fit  to  values  of  G_  as  extracted  from  the 

hn 

data  with  the  deVries  b'  fit.  Curve  1  is  the  extrapolated 
neutron-electron  interaction  slope.  Curve  2  is  the  fit  to 
the  points  indicated  by  f .  These  points  are  G„  +  , 

h  ti 

n  n 

where  G.,  is  calculated  from  the  data  using  the 
n 

Feshbach-Lomon  wave  functions  and  £G.„  are  the  relativistic 

E 

6  n 

corrections  of  Gross.  Curve  3  is  fitted  to  the  points 
shown  by  A  .  The  Partovi  wave  function  was  used  to  extract 

G_  from  these  data.  The  points  plotted  with  x  are  from 

il. 

1  ]  -3 

Drickey  and  Hand,  *  J  calculated  from  their  data  in  the  same 


way  as  those  of  curve  2. 


Figure 
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